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Dirac's method of classical analogy is employed to incorporate quantum degrees of freedom into 
modern nonequilibrium thermodynamics. The proposed formulation of dissipative quantum mechan- 
ics builds entirely upon the geometric structures implied by commutators and canonical correlations. 
A lucid formulation of a nonlinear quantum master equation follows from the thermodynamic struc- 
ture. Complex classical environments with internal structure can be handled readily. 
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Introduction. — In the context of solid state-based 
quantum information processing, external control in 
quantum optics, quantum transport through meso- and 
nanoscale structures, quantum tunneling in macroscopic 
systems, quantum Brownian motors, as well as biolog- 
ical reactions and protein folding kinetics, for example, 
quantum dissipation is an important subject. We here 
propose a general thermodynamic framework for dissipa- 
tive quantum systems in contact with classical equilib- 
rium and nonequilibrium environments. The formulation 
is guided by classical nonequilibrium thermodynamics, 
based entirely on geometric concepts and, in principle, 
applicable to arbitrarily complex quantum subsystems 
and classical environments. Important is the weak cou- 
pling of a quantum system to its environment because 
otherwise thermodynamic concepts are clearly not appli- 
cable. 

The proper arena for quantum mechanics is pro- 
vided by Hilbert spaces, that is, complex vector spaces 
equipped with inner products [HQ- ^ u the- infinite dimen- 
sional case, a Hilbert space H needs to be sufficiently rich 
to tell a nice story with all its punchlines (completeness: 
every Cauchy sequence has a limit in the space %) and 
sufficiently small to play with the tools of measure the- 
ory (separability: there exists a countable subset that is 
dense in %). Observables are densely defined self-adjoint 
linear operators on a Hilbert space T~i. As a consequence 
of self-adjointness, observables are real. The space A(H) 
of quantum observables hence is a real vector space. With 
the quantum Poisson bracket implied by the commutator 
[ , ] as a binary operation, A(TL) becomes a Lie algebra. 

For nonequilibrium systems, we are interested in the 
evolution of averages of a sufficiently rich set of observ- 
ables. It is not immediately obvious what "sufficiently 
rich" means for the characterization of the state of a 
quantum system. We here focus on the evolution of the 
density matrix p, also known as the statistical operator, 
which can be expressed in terms of an orthonormal ba- 
sis \n) of the Hilbert space % and the probabilities p n , 
namely, p — ^2 n p n \n) {n\. The density matrix character- 
izes the state of our quantum subsystem and its evolution 
determines the evolution of the averages of all observ- 
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ables. This point of view corresponds to the Schrodinger 
picture which we use throughout this letter. 

The corresponding arena for classical systems is given 
by functions or functionals on a Poisson manifold as 
observables. Classical nonequilibrium thermodynamics 
provides evolution equations for the values of observables 
in terms of Poisson and dissipative brackets as the appro- 
priate geometric structures [3|-|6|. In practice, one usually 
follows the evolution of a point i in a Poisson manifold. 
The evaluation of observables as functions or function- 
als of x is the classical counterpart of the averaging of 
quantum observables with a density matrix p. 

The classical analogy. — We here explore the possible 
mathematical setting for the formulation of nonequilib- 
rium systems with quantum degrees of freedom. The 
treatment of such systems should be inspired by the 
nonequilibrium thermodynamics of purely classical sys- 
tems. We therefore need formulations that emphasize the 
analogy between classical and quantum systems. This 
classical analogy was a major concern of P.A.M. Dirac 
who, in his pioneering textbook on The Principles of 
Quantum Mechanics devotes an entire chapter to 
the relationship between Poisson brackets and commu- 
tators which, in Chapter IV of the book |l| , are referred 
to as the quantum conditions. According to Dirac, 'We 
should thus expect to find that important concepts in 
classical mechanics correspond to important concepts in 
quantum mechanics, and, from an understanding of the 
general nature of the analogy between classical and quan- 
tum mechanics, we may hope to get laws and theorems 
in quantum mechanics appearing as simple generaliza- 
tions of well-known results in classical mechanics" (p. 84 
of Following Dirac, we here rely on his method of 
classical analogy. 

It is important to realize that, by means of the classi- 
cal analogy and simple geometric concepts, we wish to 
formulate nonequilibrium thermodynamics, not to de- 
rive it. Rigorous derivations, such as ergodic theorems, 
are almost impossible to achieve even in classical equilib- 
rium thermodynamics. Finding the proper formulation of 
nonequilibrium thermodynamics is certainly much sim- 
pler than deriving it, but no less useful. Once the proper 
thermodynamic setting is available, models for specific 
applications can be formulated readily. As in equilib- 
rium thermodynamics, closed systems, here consisting 
of quantum subsystems and their classical environments, 
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are most suited for fundamental developments and con- 
ceptual clarity. 

Reversible and irreversible structures. — A general for- 
mulation of nonequilibrium thermodynamics based on 
geometric concepts has been developed during the last 
decades . A key idea is the need for separate geomet- 
ric structures for generating reversible and irreversible 
dynamics. 

The proper structure underlying reversible dynamics is 
given known to be given by the quantum Poisson bracket 
of two observables A and B in terms of their commutator 
and Planck's constant h, 

(A,B) = ±[A,B], (1) 
in 

carefully motivated and introduced by Dirac (see p. 87 
of [H), and here written in the notation of the text- 
book p|. The rate of reversible change of the average 
(A) = tr(Ap) of a quantum observable A is given by the 
average {(A, H)) , where the observable H is the Hamil- 
tonian of the quantum system. The Poisson bracket has 
been recognized by Dirac as the common geometric struc- 
ture underlying reversible dynamics in both classical and 
quantum systems. The entire reversible structure follows 
from the noncommuting nature of quantum observables. 
Moreover, according to Heisenberg's uncertainty princi- 
ple, the averaged quantum Poisson bracket {(A,B)) oc- 
curs in the lower bound for the product of the uncer- 
tainties in two variables A and B; it certainly provides a 
most fundamental structure. 

In the same spirit, we next propose a simple geometric 
structure behind all irreversible dynamics. The essential 
structural element is ((A, Q); (B, Q)) , where Q is a suit- 
able observable that provides the dissipative coupling of 
the quantum subsystem to its environment. The bilinear 
pairing ( ; } in this structural element is the canonical 
correlation (see Eq. (4.1.12) of Q) 

(A;B) = f tY(p x Ap 1 - x B)d\ = ti-(A p B), (2) 
Jo 

where 

A p = [ P x Ap l - x d\ (3) 
Jo 

is basically the product of A and p but, as we need to 
compromise between writing p to the left or the right 
of A, we place it as a subscript. The canonical cor- 
relation is symmetric, (A;B) = (B;A) , and positive, 
(A; B) > 0. Additional convexity properties can be in- 
ferred from Lieb's theorem (see, for example, Eq. (2.120) 
of f.8j). Finally, the formula (A; 1) = (A) shows that 
canonical correlation may be considered as a pairing of 
observables and measures (density operators modifying 
the statistical operator p in the spirit of Radon- Nikodym 
derivatives), so that the dual spaces of observables and 
measures have been brought on an entirely equal footing. 



In the very same way that the reversible structure 
{(A,B)) p assists the energy to generate reversible evo- 
lution, irreversible evolution should be generated by en- 
tropy by means of the irreversible geometric structure 
((A, Q); (B, Q)) . Assuming von Neumann's logarithmic 
form of the entropy, this possibility relies on the following 
lemma, 

(In p,A p ) = (p, A). (4) 

To prove this result, we can look at arbitrary matrix ele- 
ments formed with the eigenstates of the density matrix, 

(lnp n - lnp m ) (n\ A p \m) = p m hx — / ( — J d\ 

Pm Jo \Pm J 

x (n\ A \m) = (jp n — p m ) (n\ A \m) . (5) 

By means of the lemma Q , we obtain the equally useful 
and elegant formula 

- {(A,Q);(]np,Q)) p = {(Q,(Q,A))) p . (6) 

Evolution equations. — To obtain the total energy and 
entropy of a closed system we need to know the corre- 
sponding pairs A = (A, A c ) of quantum and classical 
observables, so that one can write 

A = (A) p + A c , x . (7) 

The subscript x indicates that an observable is evaluated 
in the classical state x. Had we used probability measures 
on the state space as independent classical variables (see 
Sec. 6.3 of Q), the appearance of quantum and classical 
contributions would be even more similar, but this level 
of detail is not necessary for our purposes and most ap- 
plications. For the energy H, A is the Hamiltonian H of 
the quantum subsystem and A c is the energy H e (or E e ) 
of the environment. For the entropy 5, we choose the 
operator A as S = —fee In p, where is Boltzmann's 
constant, and A e is the entropy S c of the environment. 
We can then formulate the evolution of the average of 
any joint observable of the quantum system and its en- 
vironment in terms of the generators H and S, 

^ = P(A,H) + V(A,S), (8) 

where both the reversible Poisson contribution V and the 
irreversible dissipative contribution V consist of classical 
and quantum contributions. This formulation is the most 
natural generalization of the GENERIC framework 0- 
6] of classical nonequilibrium thermodynamics ( "general 
equation for the nonequilibrium reversible-irreversible 
coupling"). We have 

P(A,B) = {A c ,B c } x + ((A,B)) p , (9) 

where { , } is the classical Poisson bracket. Unfortu- 
nately, square and curly brackets have been used to 
distinguish between commutators and anticommutators 
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in quantum mechanics, between quantum and classi- 
cal Poisson brackets, and between Poisson and dissipa- 
tive brackets in classical nonequilibrium thermodynam- 
ics. The same symbols occur with different meanings. Al- 
though the distinction between classical Poisson brackets 
and quantum anticommutators (both expressed as curly 
brackets) as well as between classical dissipative brack- 
ets and quantum commutators (both expressed as square 
brackets) should always be clear from the context, we 
here try to facilitate the distinction by using boldface 
brackets for the classical objects. 

For the dissipative contribution to evolution we employ 
the following geometric structure, 

V(A,B) = [A ,B C ] X + [H C ,H C ]Q ((A,Q);(B,Q)) p 

- [A e ,H e ]®{(H,Q);(B,Q)} p 

- [H e ,B e ]Q((A,Q);(H,Q)) p 

+ [A e ,B e ]Q((H,Q);(H,Q)) p . (10) 

The standard classical dissipative bracket, [, ], and the 
dissipative bracket for the coupling, [, are both sym- 
metric and positive semidefinite, but they do not need 
to fulfill the usual GENERIC degeneracy expressing en- 
ergy conservation; the conservation of energy is guaran- 
teed by the form of Eq. (jTUJ) , which is actually moti- 
vated by energy conservation (that is why the last three 
compensation terms are needed). As the coupling be- 
tween the quantum system and its classical environment 
is purely irreversible, the entropy of each subsystem is 
conserved by its own reversible dynamics, {A e , S c } = 
and {(A, S)} for arbitrary A = (A, Ac). 

The dissipative structure (|TU1) can be generalized by 
summing up contributions from several different coupling 
operators Q, each of them coming with its own classical 
dissipative bracket [ , ] . The different irreversible pro- 
cesses could even be coupled. 

From our fundamental equation ([5]) , we obtain the fol- 
lowing equation for the evolution of averages in the quan- 
tum subsystem, 

^ = ((A,H)) p -[H e ,S e ]Q ((A,Q);(H,Q)) p 

+ k B [H e ,H e ]Q((Q,(Q,A))) p . (11) 

The classical environment is governed by the evolution 
equation 

^ = {A e ,H e } x -k s [A e ,H e ]Q ((Q,(Q,H))) p 

+ [A e ,S e ] x +[A e ,S e ]®((H,Q);(H,Q)) p . (12) 

Note that the average of Q in Eq. (fTTj) is not affected by 
the dissipation mediated by Q. 

Equation (|11[) for averages of the quantum subsystem 
leads us to the following quantum master equation for 



the density matrix, 

±P = -(p,H) + [H C ,S C ]Q(Q,(Q,H) P ) 

+ k B [H e ,H e ]Q(Q,(Q,p)). (13) 

The master equation (I13j) is our fundamental equation 
for open quantum systems. The existence of a master 
equation implies that the evolution of all averages is for- 
mulated in a consistent way. In view of the definition ([3]) , 
the second term in Eq. (fT3"|) will, in general, be nonlinear 
in p. This quantum nonlinearity caused by noncommu- 
tativity implies that our master equation cannot be of 
the usual Lindblad form (see, for example, Sec. 3.7 of Q 
for a discussion of nonlinear quantum master equations). 
The most natural linearization of the GENERIC mas- 
ter equation (|13[) is obtained in terms of the symmetric 
anticommutator, (Q,H) P w {(Q,H),p}/2. Note, how- 
ever, that linearizations spoil the thermodynamic struc- 
ture and are hence not recommendable. The linearized 
equation then is of the Lindblad form, where a Lindblad 
operator with real and imaginary parts proportional to 
Q and (Q,H), respectively, needs to be introduced and 
the Hamiltonian needs to be redefined. 

Instead of Eqs. (fTTj) and (fT2|) for averages, one usu- 
ally prefers to solve the evolution equations for p and x. 
The classical state variables x may be considered as spe- 
cial cases of observables, and the quantum state variable 
p can be parametrized by sufficiently many averages to 
obtain a closed set of equations. Note again that the 
counterpart of averaging quantum observables with the 
density matrix p is the evaluation of classical observables 
in a certain state x. 

Examples. — As a simple example, we consider the mo- 
tion of a particle of mass m in a potential in one dimen- 
sion. The position and momentum are given by Q and P 
with the canonical quantum Poisson bracket (Q,P) = 1, 
the potential is given by a function V(Q). We use the 
position Q as the coupling operator in Eq. (|11[) because 
friction should affect the momentum P of the particle, 
but not the position Q. The environment be a heat bath, 
the state of which be characterized by its energy H c . The 
thermodynamic properties of the bath are characterized 
by the thermodynamic relationship S C {H C ), and we in- 
troduce the bath temperature by l/T c = dS c /dH c . The 
dissipative bracket is fully determined by a friction coef- 
ficient £, 

(14) 

The equations for the first moments (Q) p , {P) p an d 
for the second moment (QQ) p agree with those for the 
widely used Caldeira-Leggett high-temperature master 
equation (see Eqs. (3.426)-(3.428) of || for ( = 2 7 m 
with a rate constant 7). In the equations for the sec- 
ond moments (PQ + QP) and (PP) however, there 
occur frictional contributions proportional to the canon- 
ical correlations (P; Q) and (P; P) rather than the 



4 



averages {PQ + QP) p /2 and (PP) p (see Eqs. (3.429)- 
(3.430) of Q)- The usual Caldeira-Leggett master equa- 
tion Q (see, for example, Eqs. (3.410) of Q) is recovered 
from the linearized GENERIC quantum master equation 
obtained after the previously mentioned approximation 
(Q, H) p w {(Q, H),p}/2 in the second term of Eq. ([I5|. 

As an example of a quantum optical master equa- 
tion, we consider the decay of a two-level system de- 
scribed in terms of the Pauli matrices ai, 03. Wc 
assume that this system is characterized by the Hamil- 
tonian H = (l/2)7kt><73, where u) is the transition fre- 
quency of the free system, and exposed to black-body 
radiation of temperature T c . The master equation for 
this system (see, for example, Eq. (3.219) of Q) is re- 
covered from the linearized GENERIC quantum mas- 
ter equation (fT3")) by choosing two coupling operators Q, 
namely <Ti and a 2 , both with equal dissipation brackets 
[H e ,H c ]® = Tchjo I (4oj) , where 70 is the spontaneous 
emission rate. To obtain the quantum master equation 
for the decay of a two-level system in its familiar form, it 
is important to note (ai, H) = —1002 and (a 2 , H) = wax 
and to introduce the rate 7 = 2~fQk-BT c /(huj). 

Quantum perspective. — In the above development, we 
have emphasized the thermodynamic perspective. Wc 
have identified the geometric structures required to for- 
mulate the thermodynamic evolution equations (j5j)- (fTU)) 
for a quantum system coupled to a classical environment, 
which are the essence of this paper. With these equations 
at hand, one can concentrate on the formulation of de- 
tailed models in terms of plausible thermodynamic build- 
ing blocks for specific real- world applications rather than 
on the theoretical derivation or justification of quantum 
master equations in simple special cases. Rather than 
heat baths only, one can treat all kinds of classical envi- 
ronments, including anisotropic and slowly changing en- 
vironments. Particularly interesting is the special case 
where the environment evolves on a time scale compara- 
ble to the time scale on which it has a significant influence 
on the quantum subsystem. 



Reaching the limit of perfectly isolated, purely re- 
versible quantum systems may not only be experimen- 
tally difficult, but also from a conceptual point of view. 
The entropy operator in Eq. ([5]) is defined if, and only if, 
p n > for all eigenstates of the density matrix. A possi- 
bly small but nonzero occupation of all quantum states 
is required for a meaningful thermodynamic description. 
Away from the idealization of a perfectly isolated system, 
the occurrence of probabilities in our theory is very nat- 
ural. As long as the limit of perfect isolation cannot be 
reached, the interpretational problems of quantum me- 
chanics seem to be absent. 

We would like to point out briefly, that the theory of 
dissipative quantum systems in Eqs. (ISf- fTU)) should be 
considered as an ideal setting for addressing a number of 
fundamental problems. As renormalization procedures 
can be related to the elimination of degrees of freedom 
on short scales, it is only natural to expect that they are 
intimately related to statistical physics. In the context 
of self-similar polymers, it has actually been shown that 
dynamical renormalization can be performed in terms 
of a renormalization of the entropy and the dissipative 
bracket [l(| • The occurrence of a black hole entropy [liT - 
IT3I ] becomes particularly interesting when it is coupled 
to other forms of entropy, including thermal entropy (see 
also [14(). The dissipative quantum mechanics (|5j)- (fTU)) 
would offer the proper framework for the description of 
such a coupling, too. 

In short, this paper suggests that there is a deep give- 
and-take relationship between quantum mechanics and 
thermodynamics. In a number of respects, quantum me- 
chanics with dissipation seems to be more natural than 
without dissipation. 

Discussions with David Taj and Jerome Flakowski mo- 
tivated me to move the topic of quantum dissipation high 
up on my priority list. I apologize for not exploring the 
pertinent literature more systematically but, as a begin- 
ner in this field, I preferred to be guided by the intrinsic 
beauty of the geometric approach. 
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